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ON THE DOUGLAS-KAZAKOV PHASE TRANSITION 

WEIGHTED POTENTIAL THEORY UNDER CONSTRAINT FOR PROBABILISTS 

THIERRY LEVY AND MYLENE MAIDA 


Abstract. We give a rigorous proof of the fact that a phase transition discovered by Douglas 
and Kazakov in 1993 in the context of two-dimensional gauge theories occurs. This phase 
transition can be formulated in terms of the Brownian bridge on the unitary group U (N) when N 
tends to infinity. We explain how it can be understood by considering the asymptotic behaviour 
of the eigenvalues of the unitary Brownian bridge, and how it can be technically approached 
by means of Fourier analysis on the unitary group. Moreover, we advertise some more or less 
classical methods for solving certain minimisation problems which play a fundamental role in 
the study of the phase transition. 
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1. Introduction 

In 1993, Douglas and Kazakov discovered that the Euclidean U(1V) Yang-Mills theory on the 
two-dimensional sphere exhibits a third order phase transition in the limit where N tends to 
infinity (see 0)- More precisely, they discovered that the free energy of this model is not a 
smooth function of the total area of the sphere. Even more specifically, they found that the 
third derivative of the free energy has a jump discontinuity when the total area of the sphere 
crosses the value n 2 . 

It turns out that the partition function of the U(./V) Yang-Mills theory on the two-dimensional 
sphere is exactly the same as that of the Brownian bridge on in a sense which we will 

explain below. In this correspondence, the area of the sphere becomes the lifetime of the bridge. 
The result of Douglas and Kazakov is thus equivalent to the existence of a phase transition 
concerning the limit as N tends to infinity of the Brownian bridge on the unitary group U(IV). 
In this limit, a certain quantity has a jump discontinuity when the lifetime of the bridge crosses 
the critical value n 2 . This is the point of view which we adopt in this paper. 

Let us give a precise statement of the main result. For all integer N > 1, let U(IV) denote 
the group of unitary N X N matrices: 

U(7V) = {U G Mjv(C) : UU* = I N }. 

Let us also denote by Hn the space of Hermitian N x N matrices: 

U N = {X € Mjv(C) : X = X*}. (1.1) 

Let us endow the real vector space FLn with the scalar product 

(X,Y) = NTr(XY). (1.2) 

Let (Xt)t> o be the corresponding linear Brownian motion on FLn- Let (Ut)t> o be the solution 
of the linear stochastic differential equation 

dU t = iU t dX t -\u t dt 

in Mjv(C), with initial condition Uq = In- The process {Ut)t >o is the Brownian motion on the 
unitary group U(1V), issued from the identity. 

For all t > 0, the distribution of Ut admits a smooth positive density with respect to the Haar 
probability measure on U(1V). We denote this density by PN,t- The main result is the following. 

Theorem 1.1. For all T > 0, the limit 

F{T) = Jim^ log Pn,t(In) (L3) 

exists. The function F : —> M is C 2 on M?j_ and C°° on M?j_ \ {vr 2 }. Moreover, the third 

derivative of F admits a jump discontinuity at the point 7r 2 . More specifically, one has 

lim F^(T) = - 7 and lim F^(T) = — 

7T 6 T—>7t 2 tt 6 

T <7T 2 T>7T 2 

In this paper, we give a complete proof of this theorem, based on the expression of the 
partition function as a series in Fourier space, that is, as a sum over Young diagrams. There is 
an alternative approach to this result, which focuses on the behaviour of the eigenvalues of the 
Brownian bridge (which we did not yet define). This point of view was taken in a previous work 
by Liechty and Wang m and it has the advantage of making it plausible, at a heuristic level, 
that there is a phase transition. We will discuss it rather informally, to the extent that it allows 
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Figure 1. The function T > F(T) (left) and its third derivative near the critical 
point 7 r 2 (right). 


one to better apprehend the nature of the phase transition. The work of Liechty and Wang, on 
the other hand, is technically arduous, and logically independent from ours. 

We also want to mention that several other third order phase transitions have been discovered 
over the last twenty years by physicists. In some recent works, Schehr and Majumdar draw a 
general picture and establish several general properties of these third order transitions. We refer 
the interested reader in particular to m and the many references therein. 

one of our main motivations in writing this 


Besides the desire to understand Theorem 1.1 


paper was to explore the concrete resolution of the minimisation problem which is at the heart 
of the study of the phase transition. In this problem, a probability measure on the real line is 
sought which minimises the sum of its energy of logarithmic self-interaction and its potential 
energy in an external field. This probability measure is moreover subject to a constraint, in our 
case that its density with respect to the Lebesgue measure should exist and not exceed 1. 

The resolution of such minimisation problems seems to share some characteristics with the 
resolution of ordinary differential equations, where more or less undisclosable manipulations are 
performed in the secrecy of one’s draft in order to identify the solution, which is then abruptly 
exhibited and proved to be the right one, using perfectly rigorous arguments which unfortunately 
give no hint as to the way in which the solution was found in the first place. 

In the present paper, we decided to disclose some of these shameful manipulations, with 
the hope that they can be fruitfully applied to other similar problems. The method which we 
explain is not at all new, we learned it in the treatise of Gakhov [ljO], and it is probably well 
known by many people who are better experts of these problems than we are. But Gakhov’s 
book is written in a rather old-fashioned language and we did not find many other accounts of 
this method: we are thus trying to write the account of it which we would have often needed 


ourselves since we are studying problems of this kind. Sections \2.4\ and \4.6\ are devoted to these 
formal derivations and we used a specific typography to emphasize that they do not have the 
same mathematical status as the rest of the paper. 

We will also try to illustrate the power of a number of tools and methods which we learned 
in the books and papers of Saff, Totik, Dragnev, Mhaskar, Kuijlaars, McLaughlin and others 
HBl El CE3 £3 IS] in order to solve these minimisation problems. In particular, we have not 
shied away from quoting several results which we found particularly useful, and likely to be 
of interest for the potential reader. To be brief, we have treated Theorem 1.1 not only as an 


interesting result in itself, but also a case study in the art of solving minimisation problems 
under constraints. 

Particle systems of the kind that will be studied in this paper are commonly referred to as 
“Coulomb gases” in the physical and mathematical literature. These systems have recently 
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received much attention, as attested for example by the works of Sandier, Serfaty and their 
coauthors (see m for a review), and also, closer to our point view, by the work of Chafai, 
Gozlan and Zitt in { 6 |, based on large deviations. In these works, very general results about the 
existence of the minimisers and their properties are stated, but only few minimisation problems 
are explicitly solvable, and even less when a constraint is added. In contrast, in this paper, we 
explore very concretely and as thoroughly as possible one very concrete example. 


Acknowledgements. We would like to acknowledge the stimulating effect which the recent 
work of Liechty and Wang m had on us. In this work, we learned among many other things 
that it is actually possible to solve rigorously the minimisation problem which is the main object 
of this paper. We also would like to thank Bernhard Beckermann for his decisive help in the 


proof of Lemma 4.5 


2. Gaussian Wigner matrices 


It will be useful for us, as a warm-up, for the sake of coherence, and in order to set up 
some notation, to review the case of Gaussian Wigner matrices in the light of the kind of 
minimisation problem which will appear in the study of the phase transition. For this now very 
classical problem, we refer mainly to the monograph [[ 2 ]. 


2.1. The Gaussian integral. For all real t > 0, the Euclidean space (%at,|(-,-)) defined 
thanks to (1.1) and (1.2) carries a Gaussian probability measure which we shall denote by 7 x,t : 

d'jNjiX) oc e“2t Tr ( x “) dx. 


Consider a function / : 1~Ln ► M which is continuous, bounded and invariant by unitary conju¬ 
gation in the sense that for all X e T~Ln and all U 6 U(IV), one has f(UXU~ 1 ) = f(X). Weyl’s 
integration formula in this context (see section 2.5.4 in |2]) asserts that 

[ f d lN,t = j, 1 ^ /v , / /(diag(Ai,..., X N ))V (Ai,..., A jV ) 2 e~i (A i + - +A ^ ) dXl , 

Jh n 11... Nl J r n 

where V(Ai,..., Ajv) = rii<i<j<iv(Ai — A?) is the Vandermonde determinant. The constant 
in front of the integral can be determined by taking / identically equal to 1 and by writing 
V(Ai,..., Ajv) as a determinant of monic Hermite polynomials. 

Introducing the empirical spectral measure 

1 N 
i=l 


the last integral can be rewritten as 


(fy 


1 1...NI 


/(diag(A 1 ,..., \ N ))e- N2jt (M d\i .^d\ N , 


%/ 2 tt 


( 2 . 1 ) 


provided we define, for every Borel probability measure /ionM with compact support, and for 
every real t > 0 , 


J* M = 


{{x,y)£R 2 :x^y} 


— log \x — y\ dfi(x)dn{y) + ^- / x 2 dy{x). 

At ./ID 


( 2 . 2 ) 


One of the keys to the intuitive interpretation of J l and the numerous similar functionals 
which we will encounter in this paper is the fact that the function (z, w) <—> — log |z — w\ is the 
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fundamental solution of the operator — on the complex plane. Accordingly, the number 
can be understood as twice the electrostatic potential energy of a unit of electric charge 
distributed according to the measure [x in the complex plane, confined in this case to the real 
line, and subject to the external potential i o 5 . One should however observe that does 
not take into account the infinite energy of self-interaction of an atom of charge with itself. 


2.2. A large deviation principle. From (2.1), one naturally expects that, when N tends to 
infinity, the eigenvalues of a typical Gaussian Hermitian matrix will be distributed so as to 
make as small as possible. More precisely, since, as it will turn out, the functional J l 

achieves its minimum at a unique diffuse probability measure, which we will denote by at, one 
expects that the probability for the empirical spectral measure fi\ to be far away from at will 
be exponentially small. This was rigorously stated and proved by Ben Arous and Guionnet in 
[ 2 ] in the framework of large deviations, as follows. 

Let us denote by A4(M) the set of Borel probability measures on M, endowed with the topology 
of weak convergence and with the corresponding Borel c-field. 

For each N > 1, let us denote by Ajv.t the distribution of the random empirical spectral 
measure jl\ of a Hermitian N x N matrix chosen under the probability measure 7jv,g 


Theorem 2 . 1 . [3j Thm 1.3] The sequence (Ajv,t)iv>i of probability measures on Ad(M) satisfies 
a large deviation principle in the scale N 2 and with good rate function 


1 


with 


TV) = 


[i i-a X’(//) + - log t- - 


-log - y\ dp(x)dp(y) + — / x 2 dy(x). 


(2.3) 


More details on the large deviation principle and on the properties of the functional X* can 
of course be found in the original paper |3j, as well as in the book [2]. 


2.3. A problem of minimisation. Let us make a couple of general remarks on the functional 
X*, most of which apply mutatis mutandis to all the functionals of the same type which we will 
consider in this paper. 

To start with, there is an issue of definition, because the first term of the right-hand side of 
is not defined for every p £ _A/f(M). This problem is solved by setting, for all (x,y) £ M 2 , 

x 2 T 

9 t(x,y) = — log |x — y\ H-——. 

The function g t is bounded below and, for all p £ Af(M), 



2 ?( m ) = / / 9t{x,y) dp(x)dp(y) 

J Jr 2 

is well defined as an element of (— 00 , 00 ]. Moreover, the map X* itself is bounded below, and 
not identically equal to 00 . Since the function gt is lower semi-continuous, the map X f itself is 
lower semi-continuous in the weak topology. A minimising sequence for X* has bounded second 
moment, so that it is tight, and X 4 attains its infimum. It is also true, although less easy to 
prove, that Tf is a strictly convex function on Ai (M) (see |18l Lemma 1.1.8] or [2( Lemma 2.6.2]), 
so that it attains its infimum at a unique probability measure. We shall denote this unique 
minimiser by at and our goal in the first part of this paper is to rediscover the shape of at- 
Before we turn to the determination of at, let us introduce some notation and collect a few 
general facts related to minimisation problems slightly more general than that of the minimisa¬ 
tion of XL 
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Consider y E A4(M) with compact support. We denote by the logarithmic potential of y, 
which is the function defined on M by 


U^x) 


log|x-y| dn{y). 


We shall also make use of the Stieltjes transform G' 1 ' of y given, for all z E C \ Supp( / u), by 


g-w = f d JM. 

Jr Z - y 


On M\ Supp(/x), one has ( U = — G^. On the support of y, and under the assumption that 
y has a sufficiently regular density, this relation takes a slightly different form, which we now 
explain. 

For any Holder continuous function ip : M —> M with compact support, we shall denote 
respectively by U^ and G^ the logarithmic potential and Stieltjes transform of the measure 
ip(x)dx. We also set 


PV / 


V’(y) 


dy = lim 


V’(y) 


x y h—^o J\y_x\>h ^ y 

the principal value of this singular integral. Then, for all x E R 


dy, 


(U^Y(x) = -Pv/ dy. 

J x-y 


(2.4) 


The following theorem describes the class of minimisation problems which we are going to 
encounter, states for each of them the existence and uniqueness of its solution, and gives a 
characterisation of this solution based on its logarithmic potential. Given a closed subset £ of 
M, we denote by A4(£) the set of Borel probability measure on £. The logarithmic energy of a 
compactly supported measure y E Af(£) is the number £(y) E (—oo,+oo] defined by 


£(id) = JJ ~ log \x - y\ dy(x)dfx{y). (2.5) 

We say that a compact subset K of C has positive capacity if there exists a probability measure 
supported by K and with finite logarithmic energy. 


Theorem 2.2 ((T8 : . Thm 1.1.3 and 1.3.3]). Let £ be a closed subset of M of positive capacity. 
Let Q : £ —> M be a lower semi-continuous function such that limM_ > . +00 (Q(a;) — log |x|) = +oo. 
Define, for any y E A4(£), 

Iq(p) = [[ -log|x-y| dy{x)dy(y) + 2 f Q(x)dy(x). (2.6) 

J Jr 2 Jr 

The infimum o/Zq over the set A4(£) is finite and it is reached at a unique probability measure 
y*, which is compactly supported. 

Moreover, assume that y is a compactly supported probability measure on £ with finite loga¬ 
rithmic energy and such that there exists a constant Fq for which 


j U^ + Q = Fq on Supp(/r), 
\ U» + Q>F q on £. 


(EL) 


Then y = y*. 


The measure y* is called the weighted equilibrium measure on the set £ in the external 
potential Q. 

The last assertion of this theorem is a powerful touchstone which allows one to check that 
a given measure is indeed the weighted equilibrium measure in a certain potential. It can be 
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understood, in electrostatic terms, by observing that + Q is the electrostatic potential jointly 
created by the distribution of charge p and the external potential Q. The fact that it is constant 
inside the support of p indicates that the charges are at equilibrium inside this support. The fact 
that it takes greater values outside this support indicates that the charges are confined within 
the support, and that taking any small amount of charge outside the support would increase 
the electrostatic energy of the system. Otherwise stated, the system (EL) is nothing but what 
physicists call the Euler-Lagrange equations associated to this problem. 

This characterisation does however not indicate how to find an expression of the measure p*. 
We will now come back to our original problem and explain how a method which we learned in 
the treatise of Gakhov |T0j , and which is also explained in Section IV. 3 of [18j, allows one to find 
a reasonable candidate for the weighted equilibrium measure on M in the potential Q(x ) = i x 2 . 
Then, in Section 2.5 we will apply Theorem 2.2 and give a rigorous proof that this candidate is 


indeed the minimiser. 

As explained in the introduction, what we explain here is not at all new, and can be regarded 
as a preparation for the much less classical minimisation problem under constraint that we will 
be facing in the second part of the paper. 


2.4. Derivation of the weighted equilibrium measure. Our goal in this section is to iden¬ 
tify ot, the unique minimiser of I 1 . As explained above, we will not care too much about rigour: 
what follows is not for the faint of heart. 

Let t > 0 be fixed. We start by making the ansatz that ut admits a density, which we denote 
by </?t, with respect to the Lebesgue measure, and that this density is regular enough, say Holder 
continuous. 

Since crt is a minimiser of It over A4(M), and assuming all the required smoothness, It does 
not vary at first order when we add to ot a small signed Borel measure of null total mass. A 
formal computation turns this observation into the assertion that the function 

i^2 / - log \x - y\Mv) dy + — x 2 (2.7) 

Jr 


must be constant on the support of at, which is nothing but the first equation of the system 
(EL) in this case. Since the quadratic external potential grows much faster than the logarithmic 
interaction potential, we expect a t to have compact support and, by symmetry, a support 
symmetric with respect to the origin. We postulatthat this support is a symmetric interval 
[—a,a]. Therefore, differentiating fl2Jl ) using \2A\ leads us to the singular integral equation 



Mv) 

x-y 


dy 


X 

— for all x E [—a, a], 


( 2 . 8 ) 


where the unknown is the function ipt. 

To solve this equation, we are going to use the Plemelj-Sokhotskyi formula (see for example 
[10] ). which asserts that if ip is Holder continuous, then for every i£l, 


Gf(x) = Pv/^M. dy^inipix), (2.9) 

J x-y 

where we use the following notation : given a function H analytic on C \ M and a point x on 
the real axis, the numbers H±(x) = lim^o H(x ± ie) are the limits from above and from below 
of H at x. 


^More details about the support will be given in Section 


2.6 
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The formula (2.9) can equivalently be written as follows: for every we have 

G%{x) + Gt(x) = 2Pv/ dy, 

J x-y 

G+(x) — Gft(x) = —2m^{x). 


( 2 . 10 ) 

( 2 . 11 ) 


Our strategy for solving (|2.8|) can be summarised as follows. 


Equation (2.11) allows us to extract g>t from the knowledge of its Stieltjes transform 
G'P*. 

Our main goal will be to find G <Pt . In view of (2.8) and (2.10), we will look for a function 
Ht analytic on C \ M which satisfies 


( H t )+{x ) + (H t )_(x) = x E [—a, a], 


( 2 . 12 ) 


and which has near infinity the behaviour of the Stieltjes transform of a probability 
measure, that is, H t (z ) = z -1 + o(z~ 2 ). This will be our candidate for G Vt . 

• Unfortunately, Equation (2.12) is not particularly easy to solve directly. As (2.11) shows, 
what is easier is to find a function analytic on C \ M with prescribed jump across the 
real axis. In order to bring (2.12) into an equation of the form (2.11), we introduce an 
auxiliary function, as we now explain. 

Let us introduce the function 

R(z) = \/z 2 — a 2 , 

analytic on C \ [—a, a], and where we choose the branch of the square root which is a positive 
real number on (a, oo). For every x E (—a, a), we have 

R±(x) = ±iV a 2 — x 2 , 

which is to say that the limits of R on the real axis from above and from below are opposite. 
This is exactly what is needed to turn the sum which appears in (2.10) into the difference which 
appears in (2.11). 


More explicitly, a function Ht satisfies (|2.12|) if and only if 

x 


H f 


R 


(H, 


(®) “ -JT (*) = ~ 2in 


\ R 


2irtV a 2 — x 2 


According to (| 2 . 11 |), a solution of this equation is given by 

12 r a 


H t z) = 


Vz 2 — a* 
2irt 


V 


'-a (z - y)\Ja 2 - y 2 


x E [—a,a]. 


dy. 


One computes 


L 


dy = 


ZTT 


~ 7T, 


and finds 


-a ( z-y)y/a 2 -y 2 \jz 2 - 

H t( z ) = -^.{z ~ \Jz 2 - a 2 ). 

The parameter a has still to be fixed. For this, we consider the behaviour of Ht near infinity: 
we have H t (z) = %z~ l + 0(z~ 2 ). On the other hand, for H t to be the Stieltjes transform of a 
probability measure, it should be equivalent to z~ 1 at infinity. This entails 

a = 2 y/i. 
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Then, for each x 6 [— 2y/t, 2 \/t], we find the function ipt, our candidate to he a solution of ( |2.8[ ) , 
by writing 

M x ) = 

This is the much expected semi-circular distribution. 

2.5. Verification of the expression of the equilibrium measure. Let us now give a proof 
of the following proposition. 


Proposition 2.3. The probability measure at given by 

d °t( x ) = —\JTt — x 2 l[_ 2 ^ 2 v^]( x ) dx - 


(2.13) 


Proof. We want to apply Theorem |2.2| For this, we need to know the logarithmic potential of 
Of. Standard computations (see for example Theorem 5.1 in [18]) show that 


is the unique minimiser over A4(M.) of the functional 27 defined by (2.3). 


U at (x) = 


1-log * 
2 

— x 2 +\x\\/x 2 —At 


J_ 

At ^ 2 


if |x| < 2 \Jt, 

— log(|x| + Vx 2 — 4t) + ^ + log2 if |x| > ‘l\Jt. 


(2.14) 


C 4f w v ^ -r 2 

Therefore, U at (x) + ^ is constant on the support of at, equal to 1 ~ 1 2 ° gt , and an elementary 
verification shows that it is larger than this constant outside the support of at- This suffices to 
guarantee that at is the minimiser of I*. □ 



Figure 2. The density of at and the graph of x i-A U at (x) + ^ for t = \. 


Since we know where X 4 achieves its minimum, we can compute this minimum, and find 


X ‘ M “ L i un{x) + it) ^ t(x) = M° gt ' 1 


In particular, the rate function appearing in Theorem 2.1 is equal to zero at at, as expected. 


2.6. The support of the weighted equilibrium measure. Let us conclude this warm up by 
mentioning that we could have approached the determination of the minimiser from a slightly 
different viewpoint, which we will adopt and illustrate in great detail in the problem of minimi¬ 
sation under constraint which constitutes the main object of this paper. 

This alternative approach relies on a series of remarkable results of Mhaskar, Saff and Totik, 
which allows one to determine before anything else the support of the weighted equilibrium 
measure. Recall from (2.5) the definition of the logarithmic energy of a compactly supported 


probability measure on C, and of the notion of positive capacity for a compact subset K of C. 


It follows from a classical result of Frostman, of which Theorem 2.2 is an elaborated version, 


that for every compact subset K of C with positive capacity, there exists a unique probability 


















10 


THIERRY LEVY AND MYLENE MAIDA 


measure u K supported by K and with minimal logarithmic energy. This probability measure is 
called the equilibrium measure of K. 

In the following theorem, and until the end of this section, we work under the assumptions 
and with the notations of Theorem 12.21 


Theorem 2.4 ([18, Thm. IV.1.5]). For every compact subset I\ of E with positive capacity, set 

MS Q (I<) = J\u“ k +Q)du K . (2.15) 

The support of the weighted equilibrium measure on E in the external potential Q is the smallest 
compact subset of E with finite logarithmic energy which minimises the functional MSq. 

The functional MSq is called the Mhaskar and Saff functional, after the names of the authors 
who introduced it. Let us emphasize that the definition which we gave for this functional is the 
opposite of the original one, simply because we find it easier to consistently minimise functionals. 

Theorem [274] alone does not suffice to identify the support of the weighted equilibrium measure, 
for there are too many compact subsets of E. However, the following theorem allows one, in 
many concrete situations, to restrict the search to a much narrower class of subsets. 


Theorem 2.5 ([15. Thm. IV. 1.10 (b)]). Consider an interval I C E such that Q is convex on 
I. Then the intersection of the support of pi* with I is an interval. 


Using Theorems 2.4 and |2.5| and an argument of symmetry, simple computations allow one 
to prove that the support of a t is [—2y/t, 2y/i\. 

One of the advantages of knowing the support of the minimising measure beforehand is the 
possibility of applying the following result. 


Theorem 2.6 ([IS, Thm. 1.3.3]). Consider a probability measure p E A4(E). Assume that v has 
finite logarithmic energy, that Supp(p) C Supp(/r*) and that U u + Q is constant on Supp(/V). 
Then v = p*. 


Concretely, if we had known that the support of the weighted equilibrium measure on M in 
the external potential x H > ^ was the interval [—2y/t,2\/i], then the first equality of (2.14) 
would have been sufficient to ensure that at was indeed the minimiser. In this situation, this 
may not seem to be a huge economy of effort, because the computation of U at is relatively easy, 
but there are more complicated situations where the computation of the potential can be rather 
challenging. 


3. Brownian motion and Brownian bridge on the unitary group 

From now on, we focus on the main object of this paper, which is to investigate the Douglas- 
Kazakov phase transition. In this section, we introduce the Brownian bridge on the unitary 
group, which is one of the models on which this transition can be studied, and we explain how, 
at a heuristic level, one can understand the origin of the Douglas-Kazakov phase transition by 
considering the behaviour of the eigenvalues of this Brownian bridge. 

3.1. The partition function of the model. Recall from the introduction that the space TLn 
of N x N Hermitian matrices is endowed with the scalar product (X,Y) = NTr(XY). The 
linear Brownian motion in 'Hat is the Gaussian process (Xt)t> o with covariance specified by the 
following relation: 

V4, B E ftjv.Vs,* > 0, E [(X a ,A)(X t ,B)] = min (s,t)(A,B). 
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The unique solution of the following linear stochastic differential equation in Mjy(C): 


dUt = iUt dXt — -Ut dt , with initial condition Uq = In, 


(3.1) 


is called the Brownian motion on the unitary group U(1V) issued from the identity. By computing 
d(UtUf), one can check that almost surely, Ut stays in U(1V) for all t > 0. 

The distribution of Ut is absolutely continuous with respect to the Haar measure on U (N) 


and we will give in (3.3) an expression of its density, which we denote by PN,t- 

Let Lf nv (\J(N)) be the space of square-integrable functions with respect to the Haar measure 
on U(1V) which are invariant by conjugation. This space admits an orthonormal basis indexed 
by the set Z^ of decreasing sequences £ = {£\ > ... > In) of IV integers, which is the Schur 
basis. For each l E Z^, the Schur function si is the quotient of the restriction to U(JV) of a 
polynomial function on Mn{ C) and an integer power of the determinant. More precisely, if U 
is a unitary matrix with pairwise distinct eigenvalues z i, • • •, zn, then 


st(U) = 


det(z i J ) i j =L ..N _ det(2 i J )jj = i...jv 


det (zf j )i,j=i...N V(zi,...,zn) 

with the notation V for the Vandermonde determinant. The value of st at matrices with multiple 
eigenvalues, for instance In, cannot be directly computed with this formula, but we have 

An) 


si{In) = Hm si{ diag(l,e* ( 

a —>-0 


a i(N— l)a^ _ V(& l ? » ■ 


an equality which is known as Weyl’s dimension formula. We can use this formula to write the 
Fourier series of the Dirac mass at In- 

V(h,...,£ N ) 


Si, = 




V(1 


S£. 


The equality above is to be understood in the following distributional sense: for every smooth 
test function / on U(V), one has /( In ) = Yleez™ V y%’ "'n') fu(N) f(U)s£(U) dU, the last integral 
being with respect to the normalised Haar measure. 

For each £ E Z^, the function is an eigenfunction of the Laplace operator A on U(JV), which 
for our present purposes can conveniently be defined as twice the generator of the Brownian 


motion defined by (3.1). The corresponding eigenvalue is the non-positive real —C 2 (£), where 


N 




i= 1 


IV- 1 


N 2 - 1 

12 


(3.2) 


For each t > 0, the Fourier series of the density of the heat kernel measure on U(!V) at time t 
(which is the distribution of Ut) can be obtained formally by applying the heat operator e? A to 
the Fourier series of the Dirac mass at the identity. As shown for example in [13J Thm 4.4 (a)], 
the resulting series is the correct one and it is normally convergent. Thus, for all U E U(JV), we 
have 




E 


We can now define the Brownian bridge as follow. Choose a real T > 0. The Brownian bridge 
of length T on the unitary group is the unique pathwise continuous stochastic process {B t ) te [o : T] 
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on U (N) such that for all n > 1, all reals 0 = < t± < ... < t n < t n+ \ = T and all continuous 

function / : \J(N) n M, 


n +1 


E[f(B t i ,..., B tn )} = ^—( f(U 1 ,...,U n )f\ pnm-U-, dU\ ... dU n , 

Z N,T JV(N) n 

with the convention Uq = U n+ \ = In, and where we have set 

„ n +1 

z n,t = / T\p N ,t i -t i -x{Ur\U i )dU 1 ...dU n . 

j U(JV)» 

This normalisation constant, or partition function, can be expressed at least in two ways: using 
the convolution property of the heat semigroup, it is equal to 


Z N,T = PN,t(In )■ 


Then, using (3.3), we find that it can also be written as 

= V r-*°2{l) V Wl, • • -Jn) 2 

v(l ,...,iV) 2 ' 


Zn.t = e 2 ' 

l&L™ 


(3.4) 


(3.5) 


The number Zn,t is our main quantity of technical interest: the Douglas-Kazakov phase 
transition expresses a failure of smothness with respect to T of its properly normalised limit as 


N tends to infinity (see Theorem 1.1). Corresponding to the two descriptions of Zn,t, there are 


two points of view on the phase transition, and two intuitive understandings of it. 

The first point of view, which was that of Douglas and Kazakov in [7], is based on the 


expression (3.5), and is focused on the behaviour as N tends to infinity of the positive measure 
on Zj of which Zn,t is the total mass. Since this is also the point of view which we adopt in 
the next sections, we will not dwell on it now. Let us simply summarise the main idea: when 
N is large, the sum is dominated by a few terms which correspond to certain vectors t which 
are close to an optimal vector, the “shape” of which depends on T and undergoes a non-smooth 
change when T crosses the critical value ir 2 . Incidentally, the reason why Douglas and Kazakov 
were interested in this quantity is that Zn.t is also the partition function of the 2-dimensional 
pure Euclidean Yang-Mills theory with structure group U(1V) on a sphere of total area T. 

The other point of view is that of the “collective field theory”, as Gross and Matytsin call it 


in DU- It is more closely related to ( |3.4[ ) , and is concerned with the behaviour of the eigenvalues 
of the Brownian bridge. Because this is not the point of view on which our technical approach 
is based, and because we find the phenomenon of transition in the behaviour of the eigenvalues 
particularly striking, we will now give a little more details about it in the two following short 
subsections. 


3.2. The eigenvalues of the Brownian motion. Before discussing the eigenvalues of the 
Brownian bridge, it is appropriate to review the more classical and easier case of the Brownian 
motion. We want to understand the global behaviour of the spectrum of a unitary matrix picked 
under the distribution of Ut- We therefore denote by 

1 N 

dU = 

2—1 

the empirical spectral measure of a unitary matrix U with eigenvalues z\,...,Zn- This is a 
probability measure on the unit circle U of the complex plane. 
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To tackle the asymptotic behaviour of the law of pu under the heat kernel measure PN,t, one 
could hope to mimic the strategy presented in Section [2] for the Hermitian case. The function 
Pnj is invariant by conjugation on U (N), so that its value at a unitary matrix U can, in principle, 
be expressed as a function of pu- Unfortunately, this is far from being as simple a function as 
the function p H > exp(—N 2 J' t (p)) which arose in the Hermitian case. 

To be clear, an analysis of the distribution of the eigenvalues of the unitary Brownian motion 
analogous to the one which we reviewed in the Hermitian case is still out of reach at the time 
of writing. A large deviation principle somewhat similar to Theorem 2.1 is very likely to hold, 
but just to identify its rate function is still an open problem. Let us nevertheless mention 
the partial results in this direction obtained, under the form of dynamical upper bounds, by 
Cabanal-Duvillard and Guionnet in [5]. 

While no large deviation principle is known, there exists a law of large numbers: the asymp¬ 
totic distribution of the eigenvalues was computed by Biane in [4j using moment techniques. To 
state his results, we use the notation tr for the normalised trace on C), the one such that 
tr (In) = 1- Given a continuous function / on HJ, we also use the notation f(U) for the result of 
the application on U of the functional calculus determined by /. Thus, tr f(U) = f v f dpu- 

Theorem 3.1 (j3j). For all t > 0, there exists a probability measure vt on HJ such that , for all 
continuous function f : HJ —> M, one has the convergence in probability 


tr f(U t ) 


TV—^oo 


/ du t 


The measure vt is determined by the fact that for all complex number z close enough to 0, 

1 du t {£) = 1 + z. 


J ml — 

This description of the measure ig is rather indirect, but using it one can determine, among 
other things, the support of For each t E [0,4], the support of v t is the interval 


e ie :\9\<l /\/— dsy 


For t > 4, the support of i ] g is the full circle U. Moreover, for every t > 0, the measure ist 
is absolutely continuous with respect to the uniform measure on HJ, and for every i / 4, this 
density is smooth and positive on the interior of the support of vt. The density of iq has a 
singularity at e m and is smooth and positive on HJ \ {e i7r }. 

For the understanding of the Douglas-Kazakov phase transition, the most important fact 
is that the support of the limiting distribution ut is equal to { 1 } when t = 0 , then grows 
symmetrically to become the full circle U at t = 4, and then stays equal to U. In other words, 
the choice which we made at the beginning for the scalar product on 'Hat (see the beginning of 
Section 3.1), and which determines the speed of the Brownian motion on U(IV), is exactly such 


that we can see, at a macroscopic scale of time, the eigenvalues of Ut progressively invade the 
whole unit circle. 

It is interesting to note that, although several quantities change qualitatively when t crosses 
the value 4, for instance the decay of the moments of the measure z^, which is polynomial for 
t < 4 and exponential for t > 4, there does not seem to be a phase transition at t = 4 in the 
same sense as there is for the Brownian bridge. 
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3.3. The eigenvalues of the Brownian bridge. Let us now turn to the behaviour of the 
eigenvalues of the Brownian bridge, the study of which is distinctly harder than in the Brownian 
case. 

It is possible to write down explicit expressions for the moments of the expected empirical 
spectral measure of the random matrix Bt , but to this day we have not been able to study their 
asymptotic behaviour directly. The analogue of Biane’s result was obtained recently by Liechty 
and Wang in the formidable paper |15j . using the point of view of determinantal processes. 
In particular, they prove the existence of a measure which plays for the Brownian bridge the 
role played for the Brownian motion by the measure i't, and they give some description of this 
measure. Note that this measure depends not only on the time t at which the Brownian bridge 
is observed, but also on the lifetime T of the bridge. Fortunately, we need not enter the details 
of Liechty and Wang’s very technical solution to understand at least at a heuristic level why the 
phase transition predicted by Douglas and Kazakov occurs. 

During the evolution of the Brownian bridge {Bt)te[o,T]i each eigenvalue starts from 1 at time 
0, wanders in the unit circle and finally comes back to 1 at time T. Collectively, the gas of 
eigenvalues undergoes an expansion on the unit circle, followed by a contraction. Since the 
processes {B t ) t £[o,T] and ( Bx-t)te[o,T] have the same distribution, we expect the evolution of 
this gas to be symmetric with respect to the time and its maximal expansion to be reached 
precisely at t = 

Now suppose that T is small. Then the time is too small for the eigenvalues to reach the 
other end of the unit circle, and the support of the limiting distribution of the eigenvalues does 
never expand enough to reach the point e J7r , which is the furthest from their starting point 1. 

If on the contrary T is large, much larger than 4, then the bridge will see its eigenvalues fill 
the unit circle almost as quickly as those of the Brownian motion (although we would expect 
the influence of the conditioning to slow this expansion slightly). Hence, at and some time 
before and some time after, the support of the limiting distribution of the eigenvalues of Bt fills 
the unit circle. 

According to this informal discussion, we expect a transition between two regimes to occur 
for a critical lifetime T of the bridge which is slightly larger than 8. The discovery of Douglas 
and Kazakov, proved by Liechty and Wang, is that the transition occurs at T = n 2 . 

We will now give another proof of the existence of the phase transition, which is closer in 
spirit to the original point of view of Douglas and Kazakov, and which is logically independent of 
Liechty and Wang’s work. Nevertheless, as mentioned in the introduction, Liechty and Wang’s 
work gave us a crucial impulsion by showing that the minimisation which is at the core of 
Douglas and Kazakov’s statement can indeed be solved rigorously. 


4. The shape of the dominant representation 

In this section, we undertake the analysis of the partition function Zn,t- We show that the 
computation of the free energy amounts to the resolution of a minimisation problem of the 
same kind as the one which we solved in the case of Wigner matrices, but with an additional 
constraint, of which we explain the nature and the origin. We then explain how to solve this 
problem, and we solve it. The computation of the free energy and the proof Theorem 0 will 
be presented in Section [5} 
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4.1. A discrete Gaussian integral. Let us start from the expression (3.5) of the partition 
function. Using (3.2), we find 


z n,t = 


e -?c 2W m 


Jn)' 




U( 1 ,...,A ) 2 


eS^ 2 - 1 ) 
V(l ,-..,A0 2 


E l (n IV — 1 \ 2 , „ „ . 

e 2 JvS i= i(h 2 ) V(£x,... ,£n) 




an expression very similar to the one which we obtained at the beginning of Section 2.1, and 
which indeed can be regarded as a discrete Gaussian integral. 

Since V(£i,... ,£n) = V(£i — ^f^, ..., £n ~the last sum is easily expressed as a function 
of the empirical measure 


W = 


We find 


Zn,t = 


~ys 

N ^ 

1=1 

e ±(N 2 -l)+N(N-l)\ogN 


J-(P — N ~ 1 ) • 
2 > 


E 


where we have set, for every Borel probability measure on 






Jt{v) = 


{(x,y)eM. 2 :x^y} 


- log \x - y\ d/j,(x)dfi(y) + — / x 2 d/x(x). 


(4.1) 


(4.2) 


Quite wonderfully, the functional Jj- is none other than Jt , as defined by (2.2). Given this 
coincidence, which of course is a reflection of a deeper correspondence between the eigenvalues 
of unitary matrices on one side and the integer vectors which index the Schur functions on the 
other side, the passage from T to is the natural consequence of our working in Fourier space. 

Recall that our main objective is to compute the free energy 

F ( T ) = Jim ^ log Z NT . 
iv—>oo 


In the scale where we are working, the contribution of the prefactor of (4.1) is easy to compute. 
Indeed, a simple computation shows that 

T '^-i)+N{N-i)]ogN\ f \ 

S ' (4.3) 


e |j(lV 2 -l)+iV(IV-l)Log]V' 

V(1 ,...,iv) 2 

We will henceforward focus on the second term, 



T 3 

“ 24 + 2 + ° 


N 


Llog •£ 


N 2 


(4.4) 




4.2. Another large deviation principle. Just as in our investigation of Wigner matrices, we 
expect the sum (4.4) to be dominated by the terms corresponding to those ts for which [if min¬ 
imises Jt- Accordingly, we expect (correctly as we will see) this phenomenon of concentration 
to give rise to a large deviation principle, and we might also expect (incorrectly as it turns out) 
that the simple relation 

F{T) = L + ^ X Q T{ai/T) = L + ^\ log T 


24 


(4.5) 
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holds, where Qt{x) = jX 2 , the functional Zq t is defined by (2.6), and <Ji/ T is the 


defined by the right-hand side of (2.13) 


mimmiser 


Let us start by stating the large deviation principle. We will then explain why (4.5) is not 
true, or at least not for all values of T. 

As always, the set A4(M) of Borel probability measures on M is endowed with the weak 
topology and the corresponding Borel (x-field. Another set of probability measures on M appears 
in the statement, namely the set 


£(M) = {p 6 A4(M) : Va, b e M, a < b, p([a, b\) < b — a} 


(4.6) 


of Borel probability measures on M which are absolutely continuous with respect to the Lebesgue 
measure, with a density not greater than 1. 


The following result can be viewed as a “discrete analogue” of Theorem 2.1 


Theorem 4.1 ( [121 Thm 2]). Let T > 0. For each N > 1, let II tv be the Borel measure on 
A4(M) defined by 

n n =Y1 ^ n2Jt{M s p 






The sequence of measures (IIjv)jv>i satisfies on A4(M) a large deviation principle in the scale 
N 2 with good rate function Zq t defined as follows: 


V/r G M(R), Z Qr {p) = 


_ ( Z Qt (h ) if /a £ £(R), 


Too 


otherwise. 


4.3. A problem of minimisation under constraint. The fact that the rate function Z Qt 
of this large deviation principle is equal to Too outside £(M) indicates that the measures fie do 
not explore the whole set A4(M). This is not surprising if we realise that the integer vectors £ 
are submitted to the constraint that their components must be pairwise distinct. Indeed, this 
constraint puts an upper bound to the concentration of mass which the measures fie can achieve, 
and we have 

f| U {fir.£eZ N l } = cm. 

p> 1 N>p 

Since £(M) is a closed subset of A4(M), the same arguments which we used after the definition 
of the functional J f (see (2.2)) imply that the functional Zq t attains its infimum on £(M) at a 
unique element of £(M), which we shall denote by /Tf. A precise statement of this fact can be 


found in Thm 2.1], and is stated below, for the convenience of the reader, as Theorem 4.4 


The existence and uniqueness of fif, combined with the large deviation principle and with 


(4.3), yields the following result, which is the key to our computation of the free energy. 


Proposition 4.2. For all T > 0, one has 

F(T) = lirn^ ^ lo S Z n,t = + ■ T e £(R)} = + (/40- 


Now, the reason why (4.5) is not true for all T is easy to explain. The problem is that oq/y, 
which is the minimiser of Zq t over the whole space A4(M), does not belong to £(M) for all T. 
In fact, a very simple computation shows that 0 \jT belongs to £(M) exactly when T < tt 2 . In 
this case, pifi = (J\it and the equality (4.5) is true, giving the correct value of the free energy. 

For T > tt 2 on the other hand, the absolute minimiser of Zq t over A4(M) does not belong to 
£(M). We are thus facing a minimisation problem of a new kind, which is the following. 
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Problem 4.3. Fix T > n 2 . Set Qt{x ) = ^x 2 . Find p*f, f/ie unique minimiser on the set £(M) 


(defined by (4.6 )) of the functional Zq t (defined by (2.6)) 


The resolution of this problem will occupy us almost until the end of this paper. Until the 
end of Section |dj we always assume that T > it 2 . 

4.4. Strategy, tools, and outline of the results. Let us start by quoting the result which 
ensures the existence and uniqueness of the constrained minimiser pf. This result is a version 
of Theorem 2.2 adapted to problems of minimisation under constraint. It is proved in a paper of 
Dragnev and Saff [9], where the authors give the first systematic study of minimisation problems 
under constraint. 

Theorem 4.4 ([9j Thm 2.1]). Let E be a closed subset o/M of positive capacity. Let Q : E —> M 
be a lower semi-continuous function such that lim| x |_ 5 . +00 (Q(x) — log|x|) = +oo. Let a be a 
Borel measure on E such that Supp(cr) = E and <r(E) > 1. Assume that the restriction of a to 
every compact set has finite logarithmic energy. 

Then the infimum of the functional Zq over the set {p 6 A4(E) : p < a} is finite and it is 
reached at a unique probability measure p *, which is compactly supported. 

Moreover, consider a compactly supported probability measure p on E such that p < a. Then 
p = p* if and only if there exists a constant Fq such that 

U ,Jj + Q < Fq p almost everywhere (a.e.) 

+ Q > Fq a — p a.e. 


(EL C 


The measure p* is called the constrained weighted equilibrium measure on E in the external 
field Q and with constraint a. 

Apart from a slight variation in the technical assumptions, the main difference between this 


result and Theorem 2.2 lies in the characterisation of the minimising measure by means of its 
potential. This new characterisation can be paraphrased by saying that 17^ + Q is equal to Fq 
on the part of the support of p where the constraint p < c r is not saturated, smaller than Fq on 
the set where the constraint p < o is saturated, and larger than Fq outside the support of p. 

In comparison with the case without constraint, the new phenomenon is that on an interval 
where p puts as much mass as it is allowed to, namely the mass given by a, it would in general 
have been energetically more efficient to put even more mass. On such an interval, U^ + Q is 
smaller than Fq and, should the constraint be released, some of the mass of p initially located 
outside this interval would migrate into it (see Figure [8] for an explicit example). 

We shall often refer to the characterisation of the constrained weighted equilibrium measure 
in terms of its potential, given by the system ( |EL e | ), as the Euler-Lagrange formulation of the 
minimisation problem. 

Let us turn to our specific problem. Considering the parity of the external potential Qt and of 
the constraint, which is the Lebesgue measure, considering the relative position of the densities 
of the absolute minimiser (Ji/t and the Lebesgue measure, and using the electrostatic intuitive 
formulation of the problem, in which the constraint can be understood as the specification of a 
finite electric capacitance of the set on which the measures live, it seems reasonable to expect 
that there will exist two reals a and f3 with 0 < a < (5 such that fif, the constrained minimiser, 
is equal to the Lebesgue measure on [—a, a], is strictly between zero and the Lebesgue measure 
on (—/?, —a) U (a, (3), and vanishes outside (—(3,f3). This is illustrated on Figure [3] below. 

According to the general principle, advertised for instance by Saff and Totik in their book 
[IS], that the first step in the determination of a minimising measure is the determination of its 
support, our first task will be to formulate this guess precisely, to prove it, and to determine 
explicitly the values of a and (3 in function of T. 
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Figure 3. A schematic representation of the measures <J\/t and hf- 


The first result which we will prove is the following. We use the notation Leb for the Lebesgue 
measure. 


Lemma 4.5. There exists two reals 0 < a < /3 such that the measure /if satisfies 
Supp(^y) = [— f3, (3] and Supp(Leb — gif) = M \ (— a, a). 


The proof of this lemma involves three ingredients. The first was suggested to us by Bernd 
Beckermann. It consists in exploiting the symmetry of the problem and transforming it into a 
problem on M + via the squaring map. The second ingredient is a duality argument, which the 
form of the statement suggests a posteriori, and which consists in studying gf and Leb — gf 
on the same footing. We learned this technique in a paper of Kuijlaars and Me Laughlin ( (13]. 
see also u Corollary 2.10]). The third ingredient is the following version of Theorem |2.5| for 
minimisation problems under constraint, which yields convexity properties of the support. 


Theorem 4.6 ( [9j Thm 2.16 (a)]). Under the assumptions of Theorem \4-4\ the intersection of 
the support of g* with any interval contained in £ and on which Q is convex is an interval. 


Let us emphasize one of the great strengths of this theorem, which is not to depend on any 
assumption on the constraining measure a. 

The actual determination of a and (3 is based on the use of the Mhaskar-Saff functional (see 
Theorem 2.4). It requires a certain amount of computation of integrals, which turn out to be 
elliptic integrals. The formulation of the result requires the definition of the complete elliptic 
integrals of the first and second kind, which are the two functions of a parameter k £ ( 0 , 1 ) 
defined by 


K(k) = 


ds 


'o \/(l ~ s 2 )(l - k 2 s 2 ) 


and E(k ) = 


1 /1 — k 2 s 2 
1 — s 2 


ds. 


(4.7) 


There is a tradition, which we follow, not to write the dependence of K{k ) and E(k) on k, and 
to write simply K and E. 


Lemma 4.7. The two reals a and j3 of which Lemma 4-5 assures the existence are uniquely 
determined, as functions of T, by the relations 

1 


P = 


2E - (1 - k 2 )K 


and a = k/3 , 


where E and K are the complete elliptic integrals defined by (4.7), and k is the unique element 
of (0,1) such that 

T = 8 EK - 4(1 - k 2 )K 2 . 

We also refer the reader to Figure [7] 
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Once a and /3 are determined, there remains to compute the restriction of pf. to the set 
J(a,/3) = [— (3, —a] U [a,/3]. For this, we use the following fact, which is a consequence of the 
result stated as Corollary 2.9 in [9]. 


Proposition 4.8. The restriction of pf to J(a,(3 ) = [—/?, — a] U [a,j3\ is (1 — 2a) times the 
weighted equilibrium mesure on the set (—oo, —a] U [a, +oo) in the external potential Qr, a = 


I+ U 


Lebn 


0 - 


At this point, the problem has been transformed into a minimisation without constraint, on 
a new subset of M, with a new potential. We will derive the solution of this problem using the 


same method which we used in Section 2.4, be it to the price of more complicated computations, 


and we will check that the solution which we obtain is indeed the minimiser using Theorem 2.6 


The last step of our proof of the phase transition consists in an explicit computation of the 
free energy of the model, and in a direct computation of the limits on the left and on the right 


at 7r" 


of the third derivative of the free energy. This will be done in Section [HJ 


4.5. The support of the weighted equilibrium potential. In this section, we prove Lemma 
14.51 and Lemma l4~7l 


Proof of Lemma f.5. In this proof, for the sake of simplicity, we shall denote the functional Zq t 
by It- 

Step 0: the minimiser is even. The fact that the potential Qt(z) = f^z 2 is even implies that 
the functional It takes the same value on a measure p and on the measure pf which is the 
image of p by the map x i —> —x. Moreover, since It is convex, this common value cannot be 
smaller than the value of It on the symmetric measure - (p + p v ). Hence, pf, the minimiser of 
It , must be a symmetric measure. 

The first idea of this proof is to exploit this symmetry and to transform the problem of 
minimisation of It into an equivalent but simpler problem on M+, via the squaring map. The 
reason why this new problem is simpler is that 0 acts as a “hard” boundary for the problem, 
and reduces the freedom with which measures are allowed to move. 

Step 1: folding up the problem. Let s : M —> M+ denote the squaring map x i-A s(x) = x 2 . 
Let us define vf = pdf o s -1 , the image measure of pf by s. We claim that uf is the unique 
minimiser of the functional 


= 


— log \u — u| du(u)du{v) + T / udu{u) 


on 


I s = ji/ G A4(M + ) : Va, b € M+,a < b,v([a,b]) < 2 (Vb-y/a) = J j . 

Indeed, a simple change of variables implies that for every p 6 £, the measure p o s~ 1 belongs 
to I s and, if p is symmetric, that the equality 

2 It{p) = If it 1 ° s ^ 1 ) 


holds. Now, every measure v in I s is the image measure by the squaring map of a unique 
symmetric measure p in £, and we have If{v) = 2 It{p) > 2 It(pt) = Ifiyf)- 

It must be noted that the squaring map plays a special role here: the logarithmic potential 
behaves particularly nicely under this transformation, and the quadratic part of the functional 
If is still the logarithmic energy. 

Let us denote by k the measure du/y/u on M + , which is the image of the Lebesgue measure 
by the squaring map. Proving the lemma is equivalent to proving that there exists two reals a 
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and b with 0 < a < b such that 


Supp(z^) = [0, b] and Supp(re — uf) = [a, +oo). 


Step 2: restricting to a compact set. We already said that our main reference for constrained 
energy problems is the paper of Dragnev and Saff [9j. In this paper, the authors work under the 
assumption that the constraint is a finite measure, an assumption which is not satisfied by the 
measure re on M_)_. Nevertheless, the proof of the existence of the minimiser (pages 243 to 246 
in 0 ) uses only the fact that the constraint is finite on compact subsets of C, which is true for 
our measure re. Moreover, as the authors explain in Remark 2.2, their proof of the existence of 
the mininriser shows that it has compact support. 

Thus, the arguments of |9] show that vtf has compact support. Let us fix M > 0, which may 
depend on T, such that the support of vtf is a subset of [0, M). Let us emphasize that, for a 
reason which shall become clear soon, we choose M large enough that it does not belong to the 
support of i/j,. Then uf is the unique minimiser of Xf on 

C S M = {v € M([0,M]) : v < K| [0 ,m]} • 


Since the measure K|[o,m] is finite, we are now exactly in the framework of the paper of Dragnev 
and Saff, and we can apply their results. 

Step 3: existence of b. Since the external potential qr{z) = 'jz to which uf is subjected is 
convex, Theorem 4.6 implies that the support of uf. is an interval. Since qx is increasing and 
the density of k is decreasing, this interval must contain 0 , for if it did not, a translation of 
vif to the left would produce another element of C s achieving a strictly smaller value of XJf. 
Altogether, this proves that the support of vf is of the form [0, b] for some b > 0, indeed b>\, 
since re([ 0 , |]) = 1 . 

Step 4: existence of a. We will now prove that the support of the difference K|[ 0 ) m] — is an 
interval of the form [ a,M ], which is to say that the measure saturates the constraint n on 
the interval [0, a). For this, we use the Euler-Lagrange formulation of the minimisation problem, 
that is, the characterisation of the minimising measure in terms of its potential (see Theorem 
4.4), and a clever duality argument which we borrow from Kuijlaars and Me Laughlin (see | [13l , 


and also [9], Corollary 2.10]). A glance at Figure [4] may help to memorise the definitions of the 
measures which we introduce. 



Figure 4. The measures re, km = ^|[o,M]; u t anc ^ v * = c( K M ~ with c = 
2y/M- 1. 
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Let us denote km = N|[o,m]- According to the Euler-Lagrange formulation of the minimisation 
problem on [0 .M] of which Uj, is the solution, there exists a constant i such that 


(EL* 


U u t + qr < i V? a.e. 
jj v t + qx > i km — v t a.e. 

Let us now set c = 2 \[M — 1 and define the probability measure 14 = \{km ~ z'jO- By merely 
rewriting the system (EL*), we see that Z 4 is the unique probability measure on [0, M] which is 
dominated by and for which there exists a constant £, the same as above, such that 

I U"* + ±(-U K ** - q T ) < *4 a.e. 

1 U v * + -q T )>-i - v* a.e. 


(EL* 


C k ' — c c 

This is exactly the Euler-Lagrange formulation of the fact that 14 is the minimiser of the 
functional 

r -U km (u)- t ~ 

— log I'll — v\ dv(u)dv{v) + 2 
l[0,M] 2 Jr+ 




V U 


2\fM — 1 


dv(u) 


on 


C s Mt = tveM([0,M]):v< 


C |[0,Af] 


2VM-1. 

This minimisation problem may look rather complicated, but all we need to prove is that the 
support of its solution 1 /* is an interval of the form [a, M], that is, an interval which contains 
M. The fact that the support of 14 contains M follows immediately from the fact that we chose 
M outside the support of v^. 

For 


In order to prove that the support of Z 4 is an interval, we apply again Theorem 4.6 


this, we need to prove that the external potential to which 14 is submitted is convex. Since qx 
is linear, it is enough to prove that U KM is concave. This is done by direct computation: one 
establishes first that for all x E [0, M], 

U KM (x ) = 2a/® log ( ] _ 2\/Mlog(M — x) + 4 Vm, 

\vM + y/x ) 

3 

then that lim^^o x? (U KM )"(x) = 0, and finally that 

d / 3 N \ \/Ma 




from which it follows that ([7 KM ) // is negative on (0, M\. This concludes the proof that JJ KM is 
concave. 

Step 5 : a 7 ^ b. We know by definition that a < b. However, we still need to rule out the 
possibility that a = b. We do this by contradiction: let us assume that a = b. Then z/ij; = N|[o, a ]) 
so that we must have a = and ^ = Leb|[_i ij. A direct computation then shows that for all 

ieR, 


( 1\, 

1 

( 1\ , 

1 

(y-2j log 

X ~ 2 

-( I+ 2j log 

* + 2 


In particular, ([//- t T) / (I) = —00 and no matter the value of T, the function U^t + Qt is strictly 
decreasing in a neighbourhood of which forbids fi,^ from being the solution of the Euler- 
Lagrange version of the minimisation problem. □ 

Let us now turn to the determination of the values of a and (3. We will use again a duality 
argument, for which we found the inspiration in [ 8 ]. 
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Proof of Lemma f.l. For all reals 0 < a < b, let us define J(a,b) = [—b,—a] U [a,b\. Let us 
denote by rf the restriction of fij, to J = J(a, j3). Let us denote by ir a the restriction to [—a, a] 
of the Lebesgue measure. By Lemma |4.5| and Proposition |4.8[ the measure rf is equal to 1 — 2a 
times the unique probability measure on = M\ (—a, a) which minimises the functional Zq t a , 
where we have set, for all 

^ _ Qt{x) + U^{x) 

Qt ’“( ) ~ 1 - 2 a 

Our goal is to determine a and /3 and this is equivalent to determining the support of rf. To do 


this, we use the functional introduced by Mhaskar and Saff and which we discussed in Section 2.6 


According to Theorem 2.4 J(a, (3) is, among all compact subsets of Sq., one which minimises the 
functional MSg Ta . From this, we deduce in particular that the function b h > MS q t a {J(a,b)) 
achieves its minimum on (a, +oo) at b = f3. An explicit computation will show that this function 
is differentiable and we find a first relation: 


d 




= 0. 


(a,ft)=(a,/3) 


(MS*) 


Since a is a boundary point of the closed set we cannot let a vary freely around a and it 
would be delicate to justify directly that J^MS Q Tc (J(a,b)) vanishes at ( a,b,c ) = (a,/3,a). To 
go around this difficulty and to get a second relation between a and /3, we use again a duality 
argument. Let us denote by 7 rg the restriction to [—/?, /3] of the Lebesgue measure, and define 
77 * = 2 /j-i ( n /3 ~ Pt) = 2 / 3—1 (Leb|j( Q| ft) — r]*). See Figure |5 for an illustration of these definitions. 





Figure 5. The measures rj* and cry*, with c = 2/3 — 1. 


From the inequalities 


it follows that 


U^t + Qt < i prp a.e. 

U^t 4- Q t > t 7T/3 — [Pf a.e. 


1 -U n P-Q t < t 

U ' 2/3-1 — 2/3-1 r l* a,e ' 

Tjn* j.- -u n P-Qr > 1 Kp _ r] „ p 

U -r 2/3-1 — 2/3-1 2/3-1 '/* d - e - 

which can conveniently be rewritten as 

ljr>* 4 . ljnp +Qt < -t n a e 

u ' 1-2/3 — 1—2/3 'I* 

tjv* 1 U np +Qr > -i JOA _ j, a P 

u ^ !_2/3 ^ 1-2/3 2/3-1 '/* 


(4.8) 


(4.9) 


(4.10) 


Since the support of 2 g-i ~ V* is the full interval [— (3, (3\, these inequalities express in Euler- 
Lagrange form the fact that 7 * is the unique minimiser of Iq t p on [—/3, (3\ . From this observation, 
it follows that the support of 77 *, which is none other than J(a,/3), achieves the minimum of 
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the Mhaskar-Saff functional MSq t/3 among all compact subsets of [—/3, /5] . In particular, taking 
differentiability for granted, 

Q Tb (J(a,b)) =0. (MS*) 

° a ’ (o,6)=(a,/3) 

Let us now compute the Mhaskar-Saff functional. In order to treat the two cases at once, let 
us choose a, b, c positive reals such that a < c <b and compute MSq t c ( J(a, &)). According to 
(2.15), we have 

(1 - 2c)MSQ T|C (J(a, b)) = (1 - 2c) J IT'M du J{afi) + f (Q + [/ Leb M) dw J(a>b) . 

It is well known (see for example Section 14 of (20j ) that the equilibrium measure of the 
compact set J(a,b) is given by 


dco 


J(a,b) 


From there one computes directly 


7r\J {b 2 — x 2 )(x 2 — a 2 ) 


1 J(a,6) (x)d'X. 


x 2 du 


+ b 2 


J(a,b) 

In order to compute the capacity of J(a, 6), one needs to compute the potential of the equilibrium 
measure uJj( a ,b)- A first step for this is to compute its Stieltjes transform 

G UJJ ( a ’ b )(z) = . ~ =, 

\J (z 2 — a 2 )(z 2 — b 2 ) 

where as usual we take the branch of the square root which is a positive real number near real 
infinity. Taking the appropriate primitive and its real part, we find that for all x 6 [—a, a], 

U Uj (o., b ) (x) = — - log f a ^ - x 2 + \J ( b 2 — x 2 )(a 2 — x 2 ) \ + - log 2. 

2 V 2 J 2 

Taking x = a, we find 

f i b 2 - a 2 

I duj( a}b) = U Uj ^ b )(a) = — - iog—-—. 

Moreover, this is also the value of U U}j ( a ’ b '< on the whole set J(a,b). Finally, we compute 

J U Leh \l-^duj J[a:b) = J L^(^) ( iLeb| hC;C] = 2 J U“ J M(x) dx. 

Since the potential of wj( a ,fe) is constant on J(a,b), hence on [a, c], we find 

J u Leh u~^dujj {ajb) = 2 U u>J< - a ’ b ) (x) dx + 2{c- a)U Uj (“’ b )(a). 

The last integral is relatively easily computed by parts, and we find, after some simplification, 

-.2 


C/ Leb |[-c,c]do; <7(ai&) = -clog 


b 2 — ar 


- 2 b K -)-E - 


Altogether, we find 


(1 - 2c)MSq t c (J(a, b)) = --log 


-7-1“ ^(° 2 + &2 ) ~ 2b [ K [ t ) “ E ( T 
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which quite remarkably does not depend on c. From there, we see in particular that the Mhaskar- 
Saff functional is indeed differentiable. Therefore, Equations (MS*) and (MS* ) are indeed satis¬ 


fied and if we denote by my (a, b ) the function on the right handside above, the numbers a and 
(3 are then solutions of the equations 

0 0 

—mr(a,/3) = 0 and —m T (a,/3) = 0. 
oa op 

Taking into account the derivatives of the complete elliptic functions, and setting k = these 
equalities entail 


k 


2/3 


i - f + 1^ 2 - rn? E = + r ' 52 ~ 2Kff+ PL? E = °- 


From these relations, one extracts 


T/3 = 4 K and P 


1 

2E-(1- k 2 )K' 


The proof is now complete. Indeed, the relation a = kp follows from our definition of k, the 
expression of P above is exactly the expected one, and the relation between T and k follows 
immediately from the two equalities above. Moreover, elementary computations allow one to 
check that the map k eA 8 EK — 4(1 — k 2 )K 2 from (0,1) to M is an increasing biiection from 
(0,1) to (tt 2 ,oo). □ 


4.6. Derivation of the constrained weighted equilibrium measure. Our problem is now 
to End the density of the measure p*, which is the restriction of the constrained minimiser p,j, 
on J(a, P) = [— P, —a] U [a, P]. According to Proposition \4.8\ rf is (1 — 2a) times the probability 
measure on J(a,P) which minimises the functional X QT, a , where for all x £ M, we have set 


„ ^ Q t (x) + U^(x) 1 ( 


— x 2 + 2a + (a 




As in Section \2.4\ we do not care much about rigour in this subsection. A rigorous proof that 
the computations we are about to do indeed yield to the minimiser will be given in the next 
subsection. 

Let us assume that rj* admits a Holder continuous density ip with respect to the Lebesgue 
measure. Let us also introduce the shorthand J = J(a,P). Recall that we denote by ir a the 
Lebesgue measure on [— a, a]. By differentiating formally the functional Zq t at pfp as we did in 
the case of Wigner matrices, or by writing the Euler-Lagrange formulation of the minimisation 
problem, we successively Gnd that the function 


x eA 


log |x 


y\ip(y) dy + 


Tx 2 

"T” 


+ H* 01 (x) 


is constant on J, and that for all x in the interior of J, 


p v fm dy = ^ + (tr-fl.) = P + bg 

Jjx-y 2 2 x + a 

In order to determine ip, we introduce the Stieltjes transform of rj*: 


G 71 (z) = 


dy*{y) 




dy. 


z-y 


j z-y 
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We are facing the same problem as in the case of Wigner matrices: we would like to know 
G 1 — G 7 L on the real axis but what we know is G T L + G^ . Let us introduce a new auxiliary 
function, which we still denote by R, and which is now analytic on C \ J: 

R(z ) = ^/(z 2 -a 2 )(z 2 -p 2 ), 

where we choose the branch of the square root which is a positive real number on (/3, oo). For 
all x 6 J, we now have 

R±(x) = ±sgn(x)ii/ (/3 2 — x 2 )(x 2 — a 2 ). 

Thus, for all x G J, 


G 71 

R 


m^(x) 

( x = - p , s T »7T- 
R+(x) 7T 


i %x + (u*°y(x) 


± ±l + ^j n R+(x) 

and, following the same reasoning as in the Wigner case, we expect the equality 


nv * ( \ * rn — 2\i 2—«2\ f Ty + ( u )(y) , 

G 1 {z) = - VO 2 -or )w -P' ) / —f - VWTW dy 

7T Jj {z-y)R+{y) 


(4.11) 


to hold. If we manage to compute the right-hand side of (|4.11|), it will be easy to compute V 7 

start with, 

( U 7ra )'(x ) = log 


using the Sokhotski formula (2.9). To start with, we have for all x G J 

x — a 


Let us introduce the function 


x + a 


r , . T , 2; — a 

LW = r + l°g—. 


which is holomorphic on C \ [—«, a], and satisfies L{x) = + {U’ Ka )\x) for all x € J. 

We are going to compute a contour integral in two different ways. Let z be a point in 
C\ [— (3, j3]. Let 7 be a contour which surrounds once both z and the interval [— /3, /3] (see Figure 
[d|). We want, to compute the integral 

f L(w ) dw 


R(w) z — w 


On one hand, we have ^7717; = O(^) near infinity, so that, by deforming 7 into a large 
circle, we see that the integral is zero. 


7 



Figure 6 . The contour 7 and its deformation. 
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On the other hand, the same contour can be deformed so as to make a small circle around z 
and to go twice along the interval [— fd, fd], once slightly below and once slightly above. We find 

rP / L+(y) + L-(y)\ dy 


Liz) 

0 = -2i7T^V4 + 


R( z ) J~/3\ R +(y) R-{y)J z-y' 

By carefully keeping track of the values of L± and R± on the intervals [—/3, —a], [—a, a] and 
[a, f3\, one finds 

L{y) dy . L(z ) . 1 dy 

— ~~ — ITT 

J ~ ( 


Jj R+{y) z ~y 

from which we conclude that 

ra 1 dy 


= —ITT 


R(z] 


—a R(y) z ~y' 


G”’ (z) = L(z) + R(z) f° A 


= log 


z — a 


z + a 


-a R(y) z ~ y 
+ - \J(z 2 - a 2 )(z 2 - fd 2 ) J 


dy 


- a 0 - y)\/(/d 2 - y 2 )(a 2 - y 2 ) 


Adding the Stieltjes transform of 7r a on one hand, setting k = ^ and performing elementary 
manipulations to the integral, we obtain our candidate for the Stieltjes transform of the full 
constrained weighted equilibrium measure ii* T : 


Ht ( z ) = ^ - irV( z2 - a 2 )( z2 - Id 2 ) 


ds 


2 ^ vv ~ yv ~ ^ o (1- ^s 2 )^/{\-s 2 ){l-k 2 s 2 ) 

This last integral is an elliptic integral of the third kind, for which one standard notation is 

ds 


(4.12) 


n {u-k) = 


’0 (1 — us 2 )y/ (1 — S 2 )(l — k 2 s 2 ) 


(4.13) 


For each k € (0,1), this is a function of v analytic in C \ [1, oo), and we can deduce from (4.12) 
our candidate for the value of the density of rf, that is of [T t , on J: 

\/x G J, ) = -^r-((H T )+ - (H t )-)(x) = V (fd 2 - x 2 )(x 2 - a 2 ) n . 

Putting everything together, we obtain the following candidate, which we denote by ipx, for the 
density of 

Tt{x) = 1 [— a,a] (x) + V(P 2 - X 2 ){x 2 - a 2 ) II 1 [-A-a]U[a lj 8](®)- ( 4 ' 14 ) 

It is interesting to note that we did not use the actual values of a and fd in our computation. 
At this point in the case of Wigner matrices, we determined the width of the support of the 
minimiser by considering the asymptotic behaviour of its Stieltjes transform. The same could 
be done here and, if we had not already rigorously determined the values of a and fd, we could 
End them by writing that Ht(z) = z -1 + 0(z~ 2 ) near infinity (see also the proof of Proposition 


5.1 below). 


The result of our derivation of is illustrated in Figure^ below. 

4.7. Verification of the expression of the equilibrium measure. We are now equipped 
with an explicit candidate for the minimiser gFp, given by (4.14) and Lemma 4.7 All we need 
to prove is that the restriction to J of the measure g>T(x)dx is, as we expect, (1 — 2a) times 
the minimiser of Iq t a on £ a . For this, according to Theorem 2.6, it suffices to check that the 
potential of tpT(x)dx is constant on J. 

The next proposition concludes the whole process of determination of the minimiser. 
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Figure 7. From left to right: k as a function of T\ a (red) and (3 (blue) as 
functions of T ; the density of the minimising measure gf as a function of T. 


Proposition 4.9. For every T > ir 2 , the unique minimiser of Iq t on £(M) is the measure 
gf. whose density pr is given by (4.14). Moreover, G^t , the Stieltjes transform of gif, is the 
function Hr given by (4.12). 


Proof. From its explicit expression (4.12), we can expand Hr near infinity. We find 

b1 £2(1 — (1 — k 2 s 2 )) 


tt ( \ T 2z/ a 2 + /3 2 

Ht ^ z) =2 z -j\!-^- 



0 -\/(l - s 2 )(l - k 2 s 2 ) 


ds I +OI 


= (|-f ) 3 + i W*® - a - *■>*))+o(^) 


and the relations between T, k and /3 imply precisely that Hr{z) = z 1 + 0(z 2 ). 

Now, the function pr defined in (4.14) is positive and smooth on M \ {±a, d=/5} and one can 


check directly that its lack of smoothness at these four points is that of square root at the origin. 
It is thus ^-Holder continuous, and the Sokhotski formula (2.9) is true at every point of M for 
its Stieltjes transform that we will denote by G^ T . Hence, G^ T and Hr have exactly the same 
jumps across the real axis. The difference between these two functions, which is analytic on C\M 
and continuous on C, is thus analytic on C, for example as a consequence of Morera’s theorem. 
By the computation above, we know that Hr is equivalent to z^ 1 near infinity, and G^ T is 
equivalent to z -1 f (pr(x)dx. Their difference is therefore entire and bounded, hence constant, 
and in fact equal to zero. This proves that Hr is the Stieltjes transform of ipr, and also that 
pr is the density of a probability measure. 

In order to see that pr is the density of the minimiser, let us denote by G^ the Stieltjes 
transform of if, which is the restriction of pr to J. The computations made in the previous 
subsection give directly that the function G^ satisfies 

V* G J, (G^+Oc) + (G^)_(x) = 2 + log . 

\ 2 x + a J 


By regularity of the function if and again by (2.9), the left-hand side is equal to 2P Vfj dy . 


Therefore, the function x i-a fj — log |x — y\ip(y)dy + —b U 7r °‘(x) is constant on J. An 
application of Proposition |4.8| concludes the proof. □ 


Figure [8] illustrates the shape of the potential of the measure gf. 
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Figure 8. The density of the minimising measure Ff and the graph of the 
function x H > U^t^x) + ?x 2 , for T slightly larger than 14. 


5. The Douglas-Kazakov phase transition 

5.1. Computation of the free energy. We now have all the tools in hand to prove our main 
result, Theorem 0 Recall from Proposition |4.2| that the free energy F, which is a function of 
T, and the regularity of which is our main concern, is the difference between an affine function 
of T and the function M defined by 

M(T) = Z t (ht)- 

Our first task is to compute M. Recall that for T > 7r 2 , the numbers k £ (0,1), j3 > 0 and 


a £ (0, jJ) are determined as functions of T by Lemma 4.7 


On the way of finding an explicit expression for the function M, that will eventually appear in 
Proposition [572], we will need to introduce the incomplete elliptic integrals of the first and second 
kind, which for each k £ (0,1) are the following analytic functions in C \ ((— oo,-l] U [1, oo])): 


F(z-k)= I 
Jo 


o y/(l - s 2 )(l - k 2 s 2 ) 


f z 11 — k 2 s 2 

ds and E(z;k)= / \ - ds, 

Jo V 1 - s 


and to find an expression of the logarithmic potential of [jJp at the points a and /3. This is what 
the next lemma does. 

Proposition 5.1. For all z £ C\ (— oo,/3], we have 

J i°g(~ - y) d^riv) = - lo s ~ a + 2 ~ 13 - 2EzF + 2KzE 0; 

- jp (i 2 ^) 2 ~ (V z 2 -a 2 + yjz 1 - ^ + l, 

(5.1) 

where log is the branch of the logarithm on C \ M_ such that to log 1 = 0. 


Proof. It suffices to check that the right-hand side of (5.1) is equal to log z + 0(z near infinity 


and that its derivative is equal to G^t. 

We now go to the main object of this subsection, namely the computation of M. 
Proposition 5.2. The function M is given by 

j i log T +\ if T < 7T 2 , 

M(r) = | 1 - ilog (^(i - fc 2 )) - |JV^(1 + * 2 ) - A/f 2 /s 2 (i - fc 2 ) 2 i/T>A. 


□ 
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The proof of this proposition is a direct computation. It can be seen as elementary, but we 
find it difficult enough to deserve a fairly detailed presentation. Moreover, we do not believe 
that this full expression has appeared in the literature before. 

Proof. Let us use the notation Q(z) = Qt(z) = jz 2 , for all z G C. To simplify the notations, 
let us denote /xf, U^t , G^t respectively by /x*, U, G. Let us also define 

9(z) = f log (z-y) d/x*(y). 

With this notation, 

M(T) = J(U + 2 Q) dy* 

and we recall that U = —whereas G' = g. 

Let us start with the expression of M in the subcritical regime T < ir 2 . In this case, we write 

M(T) = J(U + 2Q) dy* = J(U + Q) dfx* + J Q dy*. 

We know that the function U + Q is constant on the support of fx* = which is the interval 
[— ^=, ^ 7 =]. From (2.14), we know that this constant is ^(logT + 1). For the second term, a 
simple change of variable yields f Q dfx* = from where we get that M(T) = ^logT + |, 
whenever T < n 2 . 

Let us now turn to the supercritical regime T > n 2 . We will use the notation J = J(a,/3) = 
[—/5, —a] U [a, (3\. Let us split M(T ) as follows: 

/ /* pot. r pot. 

(U + 2 Q) dfx* = / (U + Q) dfx* + / Q(t) dt+ Q dfx* + / U(t) dt . 

JJ J — Ot J J —Ot 




® 


We will compute each term separately. 


1. The function U + Q is constant almost everywhere on the set J and continuous, so that it 
is constant everywhere on J. It is thus enough to evaluate it at /3. In order to evaluate U at (3, 
we use (5.1) and the relation 


U(/3) = lim — $lg(x) = — log 

x\.p 




a- 


-^ ( « 2 + / 3 2) + !- A ®. 


Thus, 


so that 


(f/ + Q)(« = -ilog^ T 2l-M(l + t 2 ) + l. 


J (U + Q) d/i* = (1 - 2a)(1/ + Q)(I3) = (2k0 - 1) P log + M {1 + k 2 ) - 1 

2 . The second term is the easiest: since Q(z) = ^z 2 and a = k/3, we have 

P Q(t) dt = \Kp 2 k z . 

J — Ot ^ 


(5.2) 

(5.3) 

(5.4) 


3. To compute the third term, we use the fact that f y 2 dfx*(y) is the coefficient of 3 i n the 
expansion of G(z) = f dfx*(y) near infinity. 












30 


THIERRY LEVY AND MYLENE MAIDA 


Recall from (4.13) the definition of the elliptic integral of the third kind. By Formula 17.7.6 


in [1], we get that for any z G \/3, oo), 


a 


n l^r-k = K + 


f3z 


y 7 ( z 2 — (3 2 )(z 2 — a 2 ) \ \ z 


(RE ( -;k) - EF ( ~-k 


We can thus rewrite (4.12) and find, for all z G \/3, oo), 


P 


G(z) = 2 EF k - 2KE k - 


P 


2 K 

T : 


1 - 


cr 




One has the following expansions in the vicinity of zero: 


F(u\ k) = u + 


1 + k 2 


6 


u 3 + 0(u 5 ) and E{u\k) = u + 


1 -k 2 
6 


« 3 + e>(u 5 ). 


Thus, 


so that 


/3 3 


G{z) = -+ £^-(l + fc 2 ) + if/3 


(1 + fc 2 ) 2 1 - fc 2 


— A; ) ) — 5 - + o ( —j 


Q4* = y(l + k 2 ) + (1 - A 2 ) 2 . 


(5.5) 


(5.6) 


4. The last term is the most delicate to compute. 

By using the symmetries of the functions U and g, and performing integrations by parts, we 
have 


/ U(t)dt = — Oft / g(t)dt = — 23? / g(t)dt = — 23? / g + (t)dt 

J—a J—a Jo Jo 

rot rot 

= 2aU(a) + 2^l tG+(t)dt = 2aU(a) + a 2 ^RG + (a) - » / f 2 G' + (f)df. 
JO Vo 

We compute t/(a) using again the fact that U + Q is constant on J: 

U(a) = (U + Q)(P) - Q(a) = log - ^(1 + 3A: 2 ) + 1. 


We then compute G+ from the (5.5): for x G [—a, a] we find 


P \ 2K ( y 7 (a 2 — x 2 )(/? 2 — x 2 ) 


+ x 


G+(x ) = 2EF. f-;k)- 2KE- I -;k ) + ^ 

\x J \ x J P \ x j 

One can check that, for x G [—a,a], 

F_ (-;k) = K- \ N ds +i f ' ds 

\ x J Ji/k \J (s 2 - l)(A’ 2 s 2 - 1) J i y 7 (s 2 - 1)(1 - & 2 s 2 ) 


P 


E_ \ —;k ) = E + 

\ x ) Jl/k 

so that 9?G+(a:) = 2Kk. 


P/ x lk 2 s 2 - 1 

s 2 — 1 


ds + i 


V k /l — fc 2 s 2 
s 2 — 1 


ds, 


Differentiating (5.7), we get that, for all x G [—a, a], 

.4 

+ 2(E - K)P 


/ \ 2iG 2 2/G 

rUG , (i) = —~x 2 H—-- - - - 

+ /3 P ^(a 2 -x 2 )(/3 2 -x 2 ) 


x 


(5.7) 


(5.8) 


y 7 (a 2 — x 2 )(/3 2 — x 2 ) 


(5.9) 
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One can check that 


[ , ^ dx = /3 (f (--k) - E (~;k)) 

Jo yj ( a 2 — x 2 )(/3 2 — x 2 ) ^ ^ « J' 


and, computing the derivative of x i-A xyj ( a 2 — x 2 )((3 2 — x 2 ), that 
„4 

-dx = 


X 


I o ^(a 2 — x 2 ){/3 2 — x 2 ) 


^/(a 2 - z 2 )(/3 2 - z 2 ) + f (a 2 + 2/3 2 )F (-;*)- ^(a 2 + /3 2 )// fc 

J Va / o \a 




Therefore, using again the relation 2 Ej3 = 1 + (1 — k 2 )K/3, we get 




K/3 


K 2 /3 2 


t 2 G' + (t)dt = -I<p 2 k 3 + -Rl + fc 2 ) + . 

3d fa 


(i-fcr-N 


There remains to put all the pieces together, and we get the expected expression for M(T). □ 

5.2. The jump of the third derivative. We finally turn to the proof of Theorem 0 

Proof of Theorem m It follows immediately from Proposition |5.2| that the function M, hence 
the free energy, is a smooth function of T on M+ \ {vr 2 }. We want to check that the function 
M^\ which is the third derivative of M, has a discontinuity of the first kind at T = 7r 2 . 

It follows from Lemma 4.7 and from the definitions (4.7) of the complete elliptic integrals 
that 

7r .. 2 


Therefore, 


lim k(T) = 0, lim K(k(T)) = K( 0) = - and lim /3(k(T)) = 

T4,7r 2 T'Itt' 2 2 t^ 7T“ 7r 


n™ M ^ = \ log ^ 2 ^ + I = 


so that M is continuous at T = n 2 . 

We will now compute the derivatives of M at 7r 2 from both sides. It turns out to be convenient 
to introduce the variable m = 1 — k 2 . Considering E and K as functions of m, one easily gets 
that, for rri 6 (0,1), 

dE 1 ._^ . OK E K 

— — = —-- (E — K) and —— =--------. 

dm 2{m — 1) dm 2m(m - 1) 2(m — 1) 

1 


The relation 2 E/3 = 1 + (1 — k 2 )I\f3 can be rewritten (5 = 

d/3 f3 


dm 4 (to — 1) 

The function M itself can be rewritten as 


(1 — mK/3 ) and 


2 E — mK 
d(Kfd) (1 — m/i/3) 2 


, so that 


dm . 


4 m(m — 1) 


M = — log/3-logm-77/3(2 — m) - m 2 K 2 (3 2 , 

2 3 12 

and we find that the derivative factorises very nicely under the form 

dM = (1 — m/i/3)(l + mKj3) . 2 _ 2m + 4) _ 

dm 24m(m — 1) 

Moreover, from the relation /3T = 4/i, one gets that 

dm m{m — l)/3 2 

Uf ~ ~ 


(1 — mK/3)(l + mK/3 )' 
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so that we finally get, that for any T > tt 2 , 

r] M R 2 

M\T) = d ^ = P -{m 2 Kp-2m + A). 

This formula is exactly the one predicted by Douglas and Kazakov when they introduced the 
model (see Equation (35) in [7]). 

From this expression of M' (T), one can check that 


lirn M'(T) 

Tin 2 


1 

2^2’ 


so that M is continuously derivable at n 2 . We keep the notation M' for the first derivative of 
M with respect to T. 

From the same formulae as above, one gets that 

dM' m/3 2 


dm 


32(m — 1) 


(1 — mK/3)( 1 + mK/3), 


so that 


This formula yields 


M"{T) 


d 2 M 

IT 2 


lirn M"(T) 

Tin 2 


so that M is twice continuously derivable at 7r 2 . 


m 2 /3 4 

32 


1 

2^’ 


Then 

= (] L h i = ( m _ 2 ( m ~ 1 ) 

dT 3 32 V1 + mK/3 (1 - mKp){l + mK/3) 

By L’Hopital’s rule, we find 

3 1 

lim M^\T) = — 77 , whereas lim M^\T) = —^, 

TU 2 K ’ 7T 6 ’ TU 2 7T 6 ’ 

so that M^ has a discontinuity of the first kind at T = ir 2 . 

Since the free energy is an affine function of T minus the function M, we find the announced 
limits on the left and on the right for the third derivatives of F. □ 
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